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Abstract — The stability of 2-k domain walls in ferromagnetic 
nanorings is investigated via calculation of the minimum energy 
path that separates a 2-k domain wall from the vortex state of a 
ferromagnetic nanoring. Trapped domains are stable when they 
exist between certain types of transverse domain walls, i.e., walls 
in which the edge defects on the same side of the magnetic strip 
have equal sign and thus repel. Here the energy barriers between 
these configurations and vortex magnetization states are obtained 
using the string method. Due to the geometry of a ring, two types 
of 2-k walls must be distinguished that differ by their overall 
topological index and exchange energy. The minimum energy 
path corresponds to the expulsion of a vortex. The energy barrier 
for annihilation of a 2k wall is compared to the activation energy 
for transitions between the two ring vortex states. 



I. Introduction 

Recent observations [£T] — fl3j in thin ferromagnetic stripes 
show magnetization configurations in which the magnetization 
makes a full 2ir turn in a localized region of the stripe, while 
the rest of the stripe is magnetized parallel to the edges of 
the stripe. These structures are stable against small applied 
external magnetic fields which makes them potentially useful 
for information storage devices 0). The same phenomenon 
has been observed in thin ferromagnetic annuli JTJ, which 
then results in the existence of a hierarchy of equally spaced 
metastable states 0. 

Reliable control of domain wall structures is crucial in the 
successful design of magnetic nanodevices 0, Q. Here we 
explore the stability of 2ir domain walls in ferromagnetic 
nanorings using the string method [8]-[ 10] . We find the energy 
barrier AE separating two metastable configurations. This 
barrier determines the escape rate from a metastable state 
through the Arrhenius law, to leading order e k B T [11]. 
The string method also gives the minimum energy path and 
transition state. 

We distinguish two types of domain walls by their winding 
number in the global (local) coordinate system to (fl). We 
compare the energy barriers that separates each type of wall 
from the ground state to the activation energy that separates 
the two vortex configurations 0, lfT2l . A current flowing 
along the axis of the ring produces a circumferential field. The 
application of this field has two effects. First, the degeneracy 
of the two vortex states is lifted. Second, the field determines 
the width of the domain wall; in its absence the 2ir walls could 
dissociate into two independent ir transverse walls. 



II. Method 

Following previous work Q, lfl2ll . we study a permalloy 
ring with the following dimensions and material properties: 
outer radius R 2 = 220 nm; inner radius Ri = 180 nm, thick- 
ness t — 2 nm, magnetization saturation M s = 8x 10 5 A/m, 
and exchange length given by ^ cx = yj ^f.^ = 5.6 nm. A 
current flowing along the axis of the ring produces a field 
H(r) = {hH c {R 1 + R 2 )/2r)e A/m. Here h = H/H c and the 
characteristic field strength at midradius is fioH c = 73.9 mT 
(for H > H c , the clockwise vortex state is no longer stable). 
The calculations were performed at h=Q. 1 . 

Precessional effects do not modify the location of the critical 
points in the energy landscape: the exponential factor in the 
Arrhenius formula is unaltered if we ignore them. We consider 
the overdamped case so that the escape trajectory follows the 
negative gradient of the energy. This is done by integrating 
only the damping term of the Landau-Lifshitz-Gilbert equation 



= -J-^-M x (M x H eff ) 

dt M s y e ' 



(1) 



Here a = 1 is the damping coefficient, 7 is the gyromagnetic 
constant, and H e ff = — Vm-EV^o is the effective magnetic 
field. The total micromagnetic energy E is the sum of the 
exchange E cx , Zeeman E% and magnetostatic terms B mag . 

The string method is necessary to calculate the minimum 
energy path between two stable states (Ma, Mb) when there 
is no a-priori knowledge of the transition state. In practice, 
the path is discretized in N+l images between Ma and 
M B denoted as M;(i) = M 4 (r,i) with i = 0,...,N. The 
images are updated using a two-step iteration procedure as 
follows: First, each image evolves using the publicly available 
micromagnetic code OOMMF [ 13 1 until the time reaches some 
interval At which we have selected to be 10 ps. This gives a 
sequence of configurations: 



M'i = M'i(r) = Mi(t) 



t+At 



dt 



dt' (2) 



Once all the M',(r) have been obtained, the second step in the 
string method is a reparametrization step used to keep these 
images equidistant. First the complete arc sn length of the 
trajectory is calculated by 



s = ) s i = s i _i + |M' i -M' !i _i| 



(3) 
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barrier for annihilation of 2n domain wall under an external 
ff] shows the configurations for images i=0,50,100. 



The arc lengths are renormalized using o! i = Si/sn- Finally 
we do a simple linear interpolation for all i along the trajectory 
so that 



Mi(t + At) = M' 



M' 



a j(»)+i a M 



M 'i(i) , i 



where j(i) is the index of the string where a'j +1 > i/N > 
a'j^y During each step we observe the magnetic energy 
Ei(t)=E„(Mi(T,t)) + E z {Mi(r,t)) + £ mag (M 4 (r, t)) as 
indicator of how far from convergence the current step is. The 
iteration process is stopped when there is no visible change in 
the function Ei(t). 

III. Annihilation of 2n domain wall. 

We now present the results of the string method to find the 
minimum energy path for destruction of a 2tt wall for the two 
types of 2tt domain wall, O = ±1. Fig. [T]and Fig. [2] present 
the string energies and configurations after relaxation of the 
string. The ft — +1, 2ir wall decays into the counterclockwise 
configuration by the expulsion of a vortex from the inner 
hole. On the other hand, the Q, = — 1, 2tt wall decays into 
the counterclockwise configuration by the expulsion of an 
antivortex. This observation shows a correspondence between 
a topological defect crossing the stripe and the signature of 
the 2tt walls being annihilated. 

For comparison purposes we provide magnitudes of the 
energy landscape of ferromagnetic nanorings obtained with 
the string method. The lowest energy barrier between the 
counterclockwise and the clockwise vortex configurations 
passes through a configuration denoted as the instanton saddle 
[12|(with O = 0); the activation energy of this event is 
3.0 x 10~ 19 J. This is consistent with our previous work on 
nanorings J5J. For a 2ir wall with f2 = 1 the decay into the 
ground state has an energy barrier equal to 1.9 x 1CP 19 J. The 
annihilation of a 2tt wall with index Q = — 1 has an energy 
barrier equal to 2.8 x 10 -19 J. This shows the energy barrier to 
annihilate a 2tt domain wall by the expulsion of a topological 
defect is comparable to that of reversal between vortex states 
by a instanton fluctuation. 




Figure 2. Segment of the ring encompassing each 2tt domain wall. Minimum 
energy path for the annihilation of 2n domain walls in consideration. (Above) 
2tt wall with topological index fi = — 1, (below) 2n wall with topological 
index Q = +1. The configurations shown correspond to the following images 
in the string (left) i=0, (center) i=50, (right) i=100. 



IV. Discussion 

For annuli with the dimensions being considered the mag- 
netization is constrained to lie in the plane of the ring; for 
the stable states the magnetization can be considered to be 
independent of the radial coordinate. The exchange energy 
E cx of a stable state is given by lfT2l : 



E R x = 



^M 2 s tf c 



27r(l + 2ft) 



2;r 



fe) de 

(5) 

where <fi(6) is the angle that the magnetization makes with the 
tangent of the ring at a given angle 8; £1 is the "winding num- 
ber" of the magnetization with respect to the local coordinate 
system. 

The difference in winding numbers (Afi) of the 2ir walls 
considered results in an exchange energy difference AEq 
between these two states. Using Eq. © the difference can 
been shown to be approximately: 



A£ n [M(r)] w 2n t x Q M^l 2 c Jn 



R2 
Ri 



An 



1.298xl0~ 19 J. 

(6) 
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The total exchange energy difference between these 2n walls 
(fi = ±1) obtained from the micromagnetic simulation results 
is 1.357 x 10 _19 J. Here we have ignored a contribution to 

this difference of terms of the form / (§f ^ dO since it is 
not a topological term. This shows a very good agreement 
between an estimate obtained from the ID model and the 
full numerical simulation. The exchange energy term is the 
biggest contributor to the difference between the total energies 
of the two domain walls: the numerical values from the 
demagnetization and Zeeman energy are 10 times smaller. The 
main point is that most of the energy difference between these 
two types of 2ir domain wall is the result of their respective 
topological windings. It is worth noting that this is a curvature 
effect as can be seen from Eq. ©: the energy difference tends 
to zero in the limit when the radii approaches infinity -i.e. in 
the limit of a straight ferromagnetic strip. 

We now consider the question of how to experimentally 
produce these two types of 2ir walls. We use the information 
contained in their global topological number u> and compare 
it to other known states of nanorings. In particular, the well- 
known "onion" state has u = 0. Since the onion is the 
remanent magnetization after saturation by an inplane uniform 
field, one can produce the = — 1 wall by applying a strong 
field in-plane followed by a circumferential field. The two 
walls will approach and form a 27r, il = — 1 wall. Changing 
the direction of either the in-plane field or the circumferential 
field will only change the final position of the 2ir wall, not its 
topological index. 

The 2tt, cj = 2 domain wall cannot be produced using only 
uniform and circumferential fields. However, we propose the 
following technique to produce that configuration in nanorings: 
apply a strong dipolar field (uj = 2) in the interior of the 
ring. This could be produced by a small current loop with its 
axis coplanar to the structure, or by bringing a magnetic tip 
close to the device. If this dipolar field is strong enough, two 
transverse walls would be produced at opposite sides of the 
ring; the magnetization vector at the centers of the wall will 
point in the opposite direction of the overall magnetization 
of the ring. Activating the circumferential field as the dipole 
strength is decreased will result in the desired configuration. 

One final question is whether resistance measurements can 
distinguish the difference between the two types of 2tt struc- 
tures reported in this work. For instance, one could attempt 
to use anisotropic magnetoresistance effect to read the overall 
winding number of the configuration. An estimate of this effect 
can be obtained by integrating J • M along a certain segment 
of the ring that spans the whole 2tt domain wall. Since J runs 
along 9 the AMR would be proportional to (Mg/M s ) which 
can be directly calculated from the 2tt domain structure. For 
two electrodes located at the top and bottom of the segment 
shown in Fig. [2] these values are 0.042 and 0.016 for the 
Q = +1 and O = — 1 walls respectively. It therefore should 
be possible to apply a current to probe the winding number 
of the 2tt domain. 



V. Conclusion 

We have presented results on the thermal annihilation of 
2ir domain walls. We differentiated between two types of 
2ir walls through their winding number in curved nanowires. 
We have observed a simple arithmetical relation between 
the topological index of the different configurations and the 
processes by which each structure decays into the ground state. 
The fact that the energy difference between the two states is 
dominated by the exchange energy allows to identify the states 
through their winding number. The transition path requires 
the motion of a singularity through the bulk: an antivortex 
destroys 51 = —1 walls, and a vortex annihilates fl = +1 
walls. Similar behavior is expected to work in linear stripes. 
Topological defects are known to play a role in certain types 
of phase transitions; here we have determined the mechanism 
by which they destroy 2tt wall structures. 

The two types of 2ir domain wall correspond to distinct 
metastable states: the greatest contribution to the energy differ- 
ence comes from the exchange energy difference. The energy 
can be directly associated with the topological signature of the 
magnetization configuration. The stability of these two states is 
comparable to the stability of the clockwise configuration. We 
suggest to use AMR noise measurements as a way to verify 
the presence of the two types of wall in magnetic nanowires 
and study their stability. Thermal fluctuations should generate 
both types of wall which could be recognized as two separate 
values of the AMR. 

Further micromagnetic exploration of the 2tt wall annihila- 
tion problem presented here can be done moving away from 
the overdamped regime by using a nonzero precessional term. 
The results presented here will be interesting to compare to 
this case. While we have presented results on a 2nm thick 
ring, we expect the observed transition states in thicker rings. 
The key parameter is the ratio of the thickness to the mean 
radius which should be less than approximately 0.1 [12 1. We 
plan to explore the energy barriers and transition states as a 
function of thickness and ring radii. 
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